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Abstract 

We model the Fock space for the electronic resonant tunneling through a double 
barrier including the coherent effects of the electron-phonon interaction. The ge- 
ometry is optimized to achieve the maximal optical phonon emission required by a 
SASER (ultrasound emitter) device. PACS numbers: 73.20.Dx, 73.40.Gk, 73.50.Rb 



The possibility of generating coherent phonons in a double barrier semicon- 
ductor heterostructure was first proposed [1] a few years ago. This is the basis 
of a SASER device [2] which transforms the electric potential energy in a single 
vibrational mode of the lattice. This is facilitated by the electronic confine- 
ment in a double barrier structure. The phonon emission appears when the 
energy of the resonant state is one quantum huJo(LO phonon energy) bellow 
the energy of the incoming electrons. As in laser devices this is enhanced if 
the first excited state of the well lies bellow the Fermi energy and becomes 
overpopulated. According to ref . [3] the emitted LO phonons decay coherently 
into a pair of LO and TA phonons the last the useful ones in a SASER device. 

In this paper we want to explore the case in which well's ground state medi- 
ates the decay of the emitter states into the collector's ones plus a phonon. 
This feature represents a resonance in the electron-phonon Fock space and is 
observed as a satellite peak in the current [4]. This resonant condition is tuned 
directly by the applied voltage and we expect that its optimization could also 
provide enough emission of primary phonons to allow for SASER operation. 
We carry out the modeling of the electronic structure and the electron-phonon 
interaction to get a minimal structure in the Fock space. Thus, the optimiza- 
tion of the phonon emission for different geometries of the device (height and 
width of the barriers, field intensity) can be discussed in simple terms. 

We consider a one-dimensional model for a double barrier including the inter- 
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Fig. 1. Each site represents a basis state in the many-body Hilbert space. The lower 
row represents localized states in different sites, the sites in black correspond to 
the barriers. The upper row represents the same sites with one phonon. The lines 
represent the nondiagonal terms in the Hamiltonian. The potential profile is also 
shown. 

action with LO phonons in the well, neglecting the effects derived from the 
accumulated charge. This will give results comparable to the 3-D case when 
Bp is small, thus limiting the number of traversal modes; or in the presence 
of a high magnetic field perpendicular to the plane of the barriers [5] which 
quantize these modes in Landau levels. We do not consider the phonon-phonon 
interaction that leads to the decay of the LO phonons. 

The Hamiltonian is a sum of an electronic contribution, a phonon contribution 
and an electron-phonon interaction term. 

Ti = Ti e + T^-p + 7~Ce-p 

Ke = Y, E j4 c 3 ~ 2 V iA4 Ck + c t c i)> 

3 j,k 

Hp = hu Hh, and H e - P = V g c k c k{bt + h). 

fc[well] A: [well] 

where c+and Cj are electron operators on site j, Ej is the site dependent 
diagonal energy and Vj fc = V5j±\^ are the hopping parameters. We assume 
that the potential drop eV is linear through the double barrier and limited to 
it. Nl and Nr are the number of sites in the left and right barriers and N w are 
those in the well, the associated lengths are Lj=Nj 2.825A. There is a single 
well state in the energy range of interest. 

Since the most important interaction between electrons and phonons in polar 
semiconductors involves longitudinal optical (LO) phonons, only one phonon 
mode with frequency Uq is considered. The electron-phonon interaction is 
limited to the well region and the coupling to the phonons is denoted with V g . 
The model is represented schematically in figure 1. 

For simplicity we restrict the problem to the case in which we have either 
or 1 phonons with no phonons in the well before the scattering process. By 
modifying V g — > V g y/n + 1 this also represents a finite temperature emission 
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Fig. 2. Inelastic current as a function of the applied voltage for different values of 
NR. 

n — > ri+1. The effective mass is taken to be 0.067 m e , the LO phonon frequency 
Tiujq = 36 meV and the value of the hopping parameter V = — 7.1018eV. 
V g is taken lOmeV which gives a typical electron-phonon interaction strength 
g = (Vg/?iu; ) 2 ~ 0.1. The barrier heights are 300meV and the Fermi energy 
e F is taken between 10 and 20 meV. 

This discrete model is solved exactly using a decimation procedure for the sites 
in the barriers and the well [6] . The leads are taken into account by adding a 
proper self-energy. The transmittances are computed from the Green's func- 
tions for the system [7]. 

Let us denote with T$~ L (T^ R ) and T^f L (T^ L ) the transmission coeffi- 
cients from left (right) to right (left) where the subscripts and 1 denote the 
number of phonons in the outgoing (first subscript) and in the incoming (sec- 
ond subscript) channel. The total current is a sum of an elastic current I e \a.nd 
an inelastic current I m (with the emission of one phonon during the scattering 
process) . These currents can be calculated from the following expressions 

I el = (2e/h) J [T R ^ L f L (e) - T^ R f R (e)]de, 
I m = (2e/h) J [T^ L f L {s) - Tfr R f R (e)]de- t 

where /l(^) and fn(e) are the Fermi functions for the left and right leads. 

For a given configuration of the system a curve of inelastic current vs. applied 
bias is obtained and its maximum value /™ ax can be extracted. Figure 2 shows 
J in -V curves as we change N R . The peaks in these curves correspond to the 
inelastic contributions to the main peak and to the satellite peak in the to- 
tal current respectively. This figure also shows that the peaks are shifted to 
higher voltages as Nr is increased. This shows a strong renormalization of the 
resonant energies due to the electrodes. In figure 3 we present i"™ 3 * vs. N R 
curves for different values of N L . These curves exhibit a maximum for J™ ax 
as a function of Nr. The optimal configurations correspond to asymmetric 
structures with wider right barriers. This can be understood by means of the 
following argument. Increasing the lifetime of the electrons in the well favors 
the electron-phonon interaction and thus increases the inelastic current. This 
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Fig. 3. Maximum inelastic current as a function of NR for Nw=20 and Ef=20 meV. 

can be done by choosing wider (or higher) barriers. In spite of this, as an effect 
of the asymmetry produced by the applied bias, the lifetime is still controlled 
mainly by the right barrier. On the other hand increasing the length of the 
barriers increases the reflectivity of the device diminishing the currents, here it 
is the left barrier which plays the main role. Then there is a trade off between 
these two effects that maximizes the phonon emission. 

In summary, we have used a simple model to show that the asymmetry in 
double barrier structures plays an important role in the ii n -V characteristics 
and to predict how it can be controlled to optimize LO phonon emission. In 
particular we show that the optimal configuration corresponds to a collector 
barrier with a length which doubles that of the emitter. 

We acknowledge finantial support from CONICET, SeCyT-U.N.C., AnPCyT 
and an international grant from Andes- Vitae-Antorchas. 
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